Spatially antisymmetric localization of matter wave in a bichromatic optical lattice 
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By direct numerical simulation of the time-dependent Gross-Pitaevskii equation using the split- 
step Fourier spectral method we study the double-humped localization of a cigar-shaped Bose- 
Einstein condensate (BEC) in a one-dimensional bichromatic quasi-periodic optical-lattice potential, 
as used in a recent experiment on the localization of a BEC [Roati et al, Nature 453, 895 (2008)]. 
Such states are spatially antisymmetric and are excited modes of Anderson localization. Where 
possible, we have compared the numerical results with a variational analysis. We also demonstrate 
the stability of the localized double-humped BEC states under small perturbation. 
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I. INTRODUCTION 

Fifty years after the prediction of Anderson localiza- 
tion |l| of electron wave in a disorder potential, the re- 
cent experimental localization [__, [3j of a non-interacting 
cigar-shaped Bose-Einstein condensate (BEC) in a quasi- 
periodic bichromatic optical-lattice (OL) 0412 and 
speckle [___| potentials has drawn much attention of re- 
search workers. The quasi-periodic bichromatic OL po- 
tential § used in the localization of a non-interacting 
BEC [|[ was formed by the superposition of two standing- 
wave polarized laser beams with incommensurate wave- 
lengths. 

Recently, there has been studies of localized BECs with 
a maximum at the center of the bichromatic OL trap, 
where the potential is a minimum pjj fl^4l^ | . A natu- 
ral extension of this phenomenon would be to investigate 
localization in more exotic situations, e.g., where a min- 
imum of the localized BEC is created at the center of 
the trap. Here, with numerical simulation of the Gross- 
Pitaevskii (CP) equation, we study the double-humped 
spatially-antisymmetric localization in an excited state 
p4H29j of a cigar-shaped BEC in a one-dimensional (ID) 
bichromatic quasi-periodic OL potential. The two humps 
appear in a single site of the OL potential. Double- 
humped structures have been created in a BEC by phase 
imprinting and other methods [30L l3l| and also have been 
studied theoretically }32M35l ]. This makes the present in- 
vestigation also of experimental interest. Anderson local- 
ization was originally predicted for the non-interacting 
system. However, in the present study on the localiza- 
tion of a double-humped BEC we also consider a weakly- 
interacting system. Although, the present localization is 
very similar to Anderson localization in a fully disordered 
potential, the bichromatic OL potential is quasi periodic 
and hence deterministic in nature. The localization con- 
sidered here is well described by the ID discrete Aubry- 
Andre model of quasi-periodic confinement [36l [37j ■ 
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FIG. 1: (Color online) (a) The bichromatic OL potential V(x) 
vs. x for different phases ao- (b) The position xo of the OL's 
central minimum vs. ao. Potential parameters in Figs. [T]to 
□ are Ai = 10, A 2 = 0.862, Ai, si = 10 and s 2 = 3. 



If the bichromatic OL potential V(x) has the sym- 
metry V(x) = V(—x), the localized states 4>{x) of the 
non-interacting BEC has the symmetry <p(x) = ±<fi(—x). 
Here we consider spatially antisymmetric localized states 
satisfying <j>{x) = —<H—x), which should be considered 
to be excited modes [2J] of Anderson localization. 

In the presence of strong disorder, the localized state 
could be quite similar to a localized state of Gaussian 
shape in an infinite potential. However, the more inter- 
esting case of localization is in the presence of a weak 
disorder when the system is localized due to the quasi- 
periodic nature of the potential [1, [23[ and not due to the 
strength of the lattice. When this happens the chemical 
potential of the system becomes comparable to the height 
of the bichromatic lattice and the localized state devel- 
ops an exponential tail. Nevertheless, the central part of 
the localized BEC is found to have a modulated Gaus- 
sian shape which allows the consideration of a variational 
approximation. 



II. ANALYTICAL CONSIDERATION 

We consider a cigar-shaped BEC under tight trans- 
verse confinement with the bichromatic OL acting along 
the axial x direction. Then it is appropriate to consider 
a ID reduction of the three-dimensional GP equation by 
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FIG. 2: (Color online) The pulse interval w of the localized 
state vs. the non-linearity g for different phase ao. 



freezing the transverse dynamics to the respective ground 
state and integrating over the transverse variables. The 
double-humped localized state of N atoms can be de- 
scribed by the dimensionless GP equation [H, HU 
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ld 2 u 
2dx^ 



g\u\ 2 u + V(x)u, 



(1) 



with normalization 



2 dx = 1, of the BEC wave 



function u = u(x,t). The spatial variable x, time t, and 
energy are expressed in transverse harmonic oscillator 
units a± = y '% / '(muS) , w" 1 and huj, where m is the mass 
of an atom and lo is the angular frequency of the trans- 
verse trap with non-linearity [38| g = 2aN/a 2 ± , and a the 
atomic scattering length. The bichromatic OL potential 
V(x) is taken as: 



V(x) = E M sin 2 (fc ; x + an), 



(2) 



with A\ = 2tt 2 si/X 2 , (I — 1,2), where A;'s are the wave- 
lengths of the OL potentials, si are their intensities, ai 
are phases, and ki = 2ix j\i the wave numbers. We take 
the phase of the first OL a\ = and that of the sec- 
ond ai = ao- In this investigation, we take the ratios 
A2/A1 = 0.862 and S2/S1 = 0.3 which arc roughly the 
same as in the experiment We further take Ai = 10, 
and si = 10. The experiment of Q employed similar 
strengths of the optical lattice to study Anderson local- 
ization in weak disorder. 

For phase difference ao = 0, potential @ has a mini- 
mum at x = 0. The position of this minimum moves to 
x = xo for ao 7^ 0. We show in Fig. [T](a) the bichro- 
matic OL potential ([2]) as ao is varied. The bichromatic 
OL potential is symmetric around x = when ao = or 
7r/2. For the phase ao in the range < ao < 7r/2, the 
OL potential is asymmetric, as shown in Fig. [T](a). In 
Fig. [T] (b) we show the position of the central minimum 
of potential © as ao is changed. 

For a variational analysis of the localized state we con- 
sider the stationary wave form 4>(x) given by u(x,t) — 



cxp(—ifit)4>(x), where /j, is the chemical potential. The 
real wave function, (f>(x), obeys the stationary equation, 

H<t>(x) + <j>"(x)/2 - g<}?(x) - V(x)<f>(x) = 0, (3) 

where the prime denotes space derivative. For ao = or 
7r/2, potential ([2]) leads to spatially-symmetric or anti- 
symmetric states confined practically to the central site 
of the quasi-periodic OL potential. When this happens, a 
variational approximation with the following Gauss-type 
ansatz is useful for the spatially-antisymmetric state 
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where w is the spatial extension of the localized state 
and will be termed pulse interval of the localized BEC 
centered at x = 0, and Af is the normalization. The 
Lagrangian of the system is given by 
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The first variational equation dL/dfi = fixes the nor- 
malization: M = 1. We use it in the following equations. 
The remaining equations dL/dw = dL/dN = yield, 
respectively, 
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(8) 



determines the pulse interval of the lo- 
The corresponding energies are given by 



Equation © 
calized state. 

E = /[((A') 2 / 2 + # 4 / 2 + Vcf) 2 ]dx = ii- 3.g/(8W27r). 

In Fig. [21 we plot the variational results for the pulse 
interval w vs. the nonlinear coefficient g according to Eq. 
([7]) . The pulse interval becomes larger as g changes from 
negative (attractive) to positive (repulsive). The reason 
is that the repulsive interaction among atoms induces a 
macroscopic repulsion between the two constituent pulses 
thus increasing the pulse interval. Compared with ao = 
0, in Fig. [2j for ao = n/2 the pulse interval is larger and 
increases faster with g. This is because, for ao = tt/2, 
the trapping is weaker relative to ao = 0. 



III. NUMERICAL RESULTS 

We perform the numerical simulation employing the 
real-time split-step Fourier spectral method with space 



FIG. 3: (Color online) Numerical (line) and variational (sym- 
bol) density \u\ 2 of the antisymmetric double-humped BEC 
vs. x for (a) g = 0, (b) g = —1, (c) g = 1, and (d) g — 2. 



FIG. 6: (Color online) Numerical density of the asymmetric 
single-humped BEC vs. x for different ao for (a) g = 1 and 
(b) g = — 1. (c) and (d) represent double-humped results for 
g = 1 and —1, respectively. 
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FIG. 4: (Color online) Numerical density \u\ 2 of localized 
states of Fig. [3] (a) vs. x in log scale. The crosses are expo- 
nential fit ~ exp( — |a;|/c0 to density, with localization length 
d = 0.17. 
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FIG. 5: (Color online) Numerical (Num) and variational (Var) 
energies and chemical potential vs. nonlinearity g. 



step 0.04, time step 0.0001. Although we use a time- 
dependent approach the localized states we calculate are 
stationary in nature. The initial input pulse is taken 
as u(x, 0) = \/2x cxp(— x 2 /2)/n 1 / 4 with a parabolic trap 
V'(x) — x 2 /2 and 5 = 0. In the beginning, during time 
evolution the parabolic trap is slowly turned off and the 
bichromatic OL is slowly turned on by increasing s\ by 
0.00001 in each time step. Successively, during time evo- 
lution we change gradually the nonlinear coefficient g by 
0.00001 in each time step to obtain the stable localized 
states. 

Figures [3] (a) , (b), (c), and (d) illustrate typical nu- 
merical and variational profiles of the localized states for 
g = 0, —1, 1 and 2 for ao = and ir/2. From Figs. [3] we 
find that the numerical densities are in good agreement 
with the variational results for small non-linearity g. If 
g is larger or the trapping is weaker, however, the lo- 
calized states develop undulating tails and occupy more 
than one OL site. In that case, density envelope of the 
localized state can not be described well by Eq. (jj) and 
the variational approximation will no longer be as good. 
We also calculated the chemical potential and energy of 
these states using Eqs. ([7]) and ([5]). The results are 
shown in Fig. [5j From Fig. [5] we find that the energy 
values are comparable to the trapping potential shown in 
Fig. Q] (a). This guarantees the interesting limit of weak 
disorder as emphasized in [H, HU . 

To see this weak disorder explicitly, we plot in Fig. 
HI the density \u\ 2 , in log scale, of the localized states 
shown in Fig. [3] (a). The long exponential tail extends 
from x = 13 to —13, whereas the central part of density 
distribution in Fig. [3] (a) contributing to normalization 
is limited between x w ±5. Wc have also shown in Fig. 
2|the exponential fit to density ~ cxp(—\x\/d) with the 
localization length [23| d — 0.17. 
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FIG. 7: (Color online) Stability dynamics \u(x, t)\ 2 vs. x and t 
for (a) the antisymmetric state with g — 1, ao — 0, [viz. Fig. 
[3] (c)] as the center of the OL trap is displaced by a small 
distance Ax = 0.1 at time t = 5, and for (b) the asymmetric 
state with g = — l,«o = 0.157T [viz. Fig. [6] (d)] as the wave 
length Ai is suddenly changed from 10 to 11 at time t = 5. 



Next we present numerical analysis for the asymmet- 
rical OL potential ([2]) with ao ^ or ^ ir/2. As shown 
in Fig. [1] (b), the position of the OL's minimum is 
no more at x = 0. On the other hand, the displace- 
ment of the OL's minimum changes periodically as ao 
increases. Because of the periodicity, we restrict our at- 
tention to < ao < 7r/2. In this case, to avoid the loss 
of accuracy in the localized state in the numerical in- 
tegration, we move the center of the initial input pulse 
and the parabolic trap to the position of the OL's min- 
imum. Thus we calculate the profiles of the localized 
states in the asymmetrical OL potential for g = 1,-1, 
and ao = 0.157T, 0.257T, and 0.347T and show the results 
in Fig. |51 In order to understand the novel effect, we 
first investigate the single-humped localized BECs with 
the asymmetrical OL potential calculated using the input 
pulse u(x,0) = exp(— x 2 /2)/n 1 / 4 . The results are shown 
in Figs. H] (a) and (b) for g = 1 and —1, respectively. 
It can be seen that the atom density envelopes of the 
single-hump localized BEC arc asymmetrical. Compared 
with Fig. [5] (a), Fig. [5] (b) presents narrower and larger 
atom densities because of the focusing nature of the BEC 
interaction (g = —1). When ao = 0.347T, the pulse width 
of the localized state is larger compared to the widths 
for ao = 0.157T, or 0.257T. The reason is that the overall 
trapping is weaker as ao increases. 

The same results are next obtained for the double- 
humped localized BECs which are shown in Figs. [5] (c) 
and (d) for g = 1 and g = — 1, respectively, using the ini- 
tial input pulse u(x, 0) = \/2xcxp(~ x 2 /2)/ir 1 / i . These 
two Figs, show the asymmetry of the two humps for 
the double-humped localized BECs. The right hump of 
the localized states is higher than its left hump because 
trapping on the left is stronger. When ao = 0.347T, the 
difference in the height of the two humps and the pulse 
interval of the localized state are larger compared with 
those of ao = 0.25-7T or 0.157T. These results arc in agree- 
ment with single-humped localized BECs illustrated in 
Figs. M (a) and (b). 

One of the most important issues is the stability of the 
double- humped localization against perturbation. First, 
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FIG. 8: (Color online) The numerical (solid line) and vari- 
ational (dotted line with symbol) density \u\ 2 vs. i of the 
double-humped BEC for Ai = 6 and 8, and g = 0, si = 
10, S2 = 3, A2 = 0.862Ai, and ao = 0. The corresponding 
bichromatic OL potential is also shown in arbitrary units. 



we investigate the stability of the antisymmetrical lo- 
calization. In order to examine whether the predicted 
localized state is stable, we introduce a small perturba- 
tion from equilibrium point by displacing the center of 
the OL trap by a small distance Ax = 0.1 at t = 5. 
The numerical simulations are shown in Fig. [JJ (a) for 
g = l,ao = 0. The double- humped localization oscillate 
around the new center of the OL trap, and the symmetry 
of the two humps is broken after the small perturbation is 
introduced. The double-humped localization exists dur- 
ing a large time interval for small perturbation. We also 
checked the stability of the states for other values (nega- 
tive) of g and for other types of perturbation by changing 
the value of the wave length Ai and the antisymmetrical 
localized state was found to be stable. Next we stud- 
ied the stability of the asymmetric localized states. The 
stability of one such state is shown in Fig. [7J (b) for 
ao = 0.157T and g — —1 when Ai was suddenly changed 
from 10 to 11 at time t = 5. The double-humped state is 
again found to be stable against small perturbation. 

So far we considered a fixed OL wave length: Ai = 10. 
To get insight into the effects of the OL on the station- 
ary localized states, we now investigate the properties 
of the system when the wave lengths A; arc smaller. In 
realistic experiment, the dimensionless wavelength can 
be adjusted by the harmonic trap because Ai is related 
to the transverse harmonic-oscillator length. The wave 
length Ai controls the dimensionless heights of potential 
@, so it is meaningful to investigate the effects of Ai 
on the stationary localized states. In this part, the same 
parameters of the numerical integration are selected as 
those in the preceding calculation, (viz. si = 10, S2 = 3, 
A 2 = 0.862Ai, and a = 0). In Fig. [5]we plot |m| 2 vs. x 
for g = and different Ai. A smaller Ai leads to a local- 
ized state with a smaller pulse interval w. We have also 
compared the numerical results with variational analysis 
for small Ai. The variational results are shown by the 
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chain of symbols in Fig. [8j 

IV. SUMMARY 

Using the numerical and variational solution of the 
GP equation, wc studied the stationary localization of 
a double-humped cigar-shaped BEC in a bichromatic 
quasi-periodic ID OL potential. The bichromatic OL po- 
tential is generated by superposing two OL potentials in 
the form of sine waves. Such a bichromatic OL potential 
is symmetric around the center at x — 0, consequently, 
the density of the localized state also possesses the same 
symmetry. In the presence of a phase difference between 
the two OL components the above symmetry is broken 



and we analyze this symmetry breaking in case of single- 
humped and double-humped states. Here we also study 
the effect of a weak atomic interaction (both attractive 
and repulsive) on the profile of the localized states. The 
localized double-humped states were found to be dynam- 
ically stable under small perturbations. We hope that 
the present work will motivate new studies, specially ex- 
perimental ones on Anderson localization in the form of 
double-humped states. 
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